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^ . Abstract 

oo ; 

\Q , Observational evidence along with theoretical arguments which call for non-baryonic dark mat- 
O ' ter are reviewed. A brief summary of the dark matter session is included. 

O 

Q\ . There is increasing evidence that relative to the visible matter in the Universe, which is 

^ . in the form of baryons, there is considerably more matter in the Universe that we don't see 
Here, I will review some of the motivations for dark matter in the Universe. The best 
2 ! observational evidence is found on the scale of galactic halos and comes from the observed flat 
^ ! rotation curves of galaxies. There is also mounting evidence for dark matter in elliptical galaxies 
as well as clusters of galaxies coming from X-ray observations of these objects. Also, direct 
evidence has been obtained through the study of gravitational lenses. In theory, we believe there 
^ I is much more matter because 1) inflation tells us so (and there is at present no good alternative 
I to inflation) and 2) our current understanding of galaxy formation only makes sense if there is 
more matter than we see. One can also make a strong case for the existence of non-baryonic 
dark matter in particular. The recurrent problem with baryonic dark matter is that not only 
is it very difficult to hide baryons, but the standard model of primordial nucleosynthesis would 
have to discarded if all of the dark matter is baryonic. Fortunately, as will be covered at length 
in these proceedings, there are several attractive alternatives to baryonic dark matter. 

Before embarking on the subject of dark matter, it will be useful to review the relevant 
quantities from the standard big bang model. In a Friedmann-Robertson- Walker Universe, the 
expansion rate of the Universe (the Hubble parameter) is related to the energy density p and 
curvature constant k by 

. fR\^ SttG k 
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assuming no cosmological constant, where k = ±1, for a closed, open or spatially flat Universe, 
and R is the cosmological scale factor. When k = 0, the energy density takes its "critical" value, 
p = Pc = 3H^/8nG = 1.88 x lO'^^hl g cm'^ where K = Ho/WO km s"^ Mpc"^ is the scaled 
present value of the Hubble parameter. The cosmological density parameter is defined by 
Q = p/pc and by rewriting eq. (0) we can relate k to Q and H by 

= (fi - (2) 

so that k = +1, —1, corresponds to^7>l,<l,= l. 

In very broad terms, observational limits on the cosmological parameters are: 0.2 < Q < 2 
and 0.4 < /iq < 1.0 [@. The cosmological density is however sensitive to the particular scale 
being observed (at least on small scales). Accountably visible matter contributes in total only 
a small fraction to the overall density, giving fly ~ -003 — .01. In the bright central parts of 
galaxies, the density is larger Q ~ 0.02 — 0.1. On larger scales, that of binaries and small groups 
of galaxies, Q ~ 0.05 — 0.3. On even larger scales the density may be large enough to support 
Q ~ 1.0. Though there are no astronomical observations to support Q > 1, limits based on the 
deceleration of the Universe only indicate that Q <2. 

The age of the Universe is also very sensitive to these parameters. Again, in the absence 
of a cosmological constant we have, 

Hotu= f\i-n + n/x)-^/^dx (3) 

Jo 

For tu > 13 Gyr, nhl < 0.25 if K > 0.5 and nhl < 0.45 if K > 0.4. While for tu > 10 Gyr, 
Qhl < 0.8 if K > 0.5 and Qh^ < 1.1 if K > 0.4. 

There is, in fact, good evidence for dark matter on the scale of galaxies (and their halos). 
Assuming that galaxies are in virial equilibrium, one expects that by Newton's Laws one can 
relate the mass at a given distance r, from the center of a galaxy to its rotational velocity 

M(r) oc v^v/Gn (4) 

The rotational velocity, v, is measured 0, ^, ^ by observing 21 cm emission lines in HI regions 
(neutral hydrogen) beyond the point where most of the light in the galaxy ceases. A compilation 
of nearly 1000 rotation curves of spiral galaxies have been plotted in as a function of r for 
varying brightnesses. If the bulk of the mass is associated with light, then beyond the point 
where most of the light stops M would be constant and oc 1/r. This is not the case, as 
the rotation curves appear to be fiat, i.e., v ~ constant outside the core of the galaxy. This 
implies that M oc r beyond the point where the light stops. This is one of the strongest pieces 
of evidence for the existence of dark matter. Velocity measurements indicate dark matter in 
elhptical galaxies as well 0. 

Galactic rotation curves are not the only observational indication for the existence of 
dark matter. X-ray emitting hot gas in elliptical galaxies also provides an important piece 
of evidence for dark matter. As an example, consider the large elliptical M87. The detailed 
profiles of the temperature and density of the hot X-ray emitting gas have been mapped out 
1^. By assuming hydrostatic equilibrium, these measurements allow one to determine the 
overall mass distribution in the galaxy necessary to bind the hot gas. Based on an isothermal 
model with temperature kT = 3keV (which leads to a conservative estimate of the total mass). 
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Fabricant and Gorenstein predicted that the total mass out to a radial distance of 392 Kpc, 
is 5.7 X IO^^Mq whereas the mass in the hot gas is only 2.8 x IO^^Mq or only 5% of the total. 
The visible mass is expected to contribute only 1% of the total. The inferred value of Q based 
on M87 would be ~ 0.2. 

M87 is not the only example of an elliptical galaxy in which X-ray emitting hot gas is 
observed to indicate the presence of dark matter. At this meeting, Foreman showed several 
examples of ellipiticals with large mass to light ratios. For example in the case of N4472, while 
the optical observations go out to 25 kpc, the X-ray gas is seen out to 75 kpc, indicating M/L's 
of about 60 at 70 kpc and up to 90 at 100 kpc. Similar inferences regarding the existence of 
dark matter can be made from the X-ray emission from small groups of galaxies [|10|, |lT|. 

On very large scales, it is possible to get an estimate of fl from the distribution of peculiar 
velocities. On scales. A, where perturbations, S, are still small, peculiar velocities can be 
expressed as f ~ HX6Q^'^. On these scales, measurements of the peculiar velocity field 



from the IRAS galaxy catalogue indicate that indeed fl is close to unity [|TH|. Another piece 



of evidence on large scales, is available from gravitational lensing [|14|. The systematic lensing 
of the roughly 150,000 galaxies per deg^ at redshifts between 2; = 1 — 3 into arcs and arclets 
allow one to trace the matter distribution in a foreground cluster. Van Waerbeke discussed 
recent results of weak gravitational lensing looking at systems, which if virialized, have mass 
to light ratios in the range 400-1000 and correspond to values of fl between 0.25 and 0.6 [p!5[] . 
Unfortunately, none of these observations reveal the identity of the dark matter. 

Theoretically, there is no lack of support for the dark matter hypothesis. The standard 



big bang model including inflation almost requires that ^1 = 1 |T^. The simple and unfortunate 
fact that at present we do not even know whether Q is larger or smaller than one indicates that 
we do not know the sign of the curvature term, further implying that it is subdominant in Eq. 

(13) 

k 8nG 

< P (5) 

In an adiabatically expanding Universe, R ~ where T is the temperature of the thermal 
photon background. Therefore the quantity 

k = —— < — - < 2 X 10"^^ 6 
R^T^ 3T2 ^ ^ 

is dimensionless and constant in the standard model. The smallness of k is known as the 
curvature problem and can be resolved by a period of inflation. Before inflation, let us write 
R = Ri, T = Ti and R ~ T~^. During inflation, R ~ ~ e^*, where H is constant. After 
inflation, R = Rj ^ Ri but T = Tj = < where Tr is the temperature to which the 
Universe reheats. Thus R 7^ T^^ and k —* is not constant. But from Eqs. (Q) and if 
A; — s> then Q ^ 1, and since typical inflationary models contain much more expansion than is 
necessary, Q becomes exponentially close to one. 

If this is the case and f2 = 1, then we know two things: Dark matter exists, since we don't 
see = 1 in luminous objects, and most (at least 90%) of the dark matter is not baryonic. 
The latter conclusion is a result from big bang nucleosynthesis 1171 , which constrains the 
baryon-to-photon ratio rj = hb /n-^ to 

1.4 X 10"^° < < 3.8 X 10"^° (7) 



3 



which corresponds to a hmit on fig 



0.005 <Qb< 0.09 (8) 

for 0.4 < ho < 1.0. Thus 1 — VLb is not only dark but also non-baryonic. I will return to big 
bang nucleosynthesis below. 

Another important piece of theoretical evidence for dark matter comes from the simple 
fact that we are here living in a galaxy. The type of perturbations produced by inflation [|19| 
are, in most models, adiabatic perturbations {5p/ p oc 6T/T), and I will restrict my attention 
to these. Indeed, the perturbations produced by inflation also have the very nearly scale-free 
spectrum described by Harrison and Zeldovich [^. When produced, scale-free perturbations 
fall off as ^ oc /^^ (increase as the square of the wave number). At early times 6p/p grows 
as t until the time when the horizon scale (which is proportional to the age of the Universe) 
is comparable to Z. At later times, the growth halts (the mass contained within the volume 
has become smaller than the Jean's mass) and ^ = S (roughly) independent of the scale I. 
When the Universe becomes matter dominated, the Jean's mass drops dramatically and growth 
continues as — oc i? ~ 1/T. For an overview of the evolution of density perturbations and the 



resulting observable spectrum see The transition to matter dominance is determined by 
setting the energy densities in radiation (photons and any massless neutrinos) equal to the 
energy density in matter (baryons and any dark matter). For three massless neutrinos and 
baryons (no dark matter), matter dominance begins at 

= 0.22mBV (9) 

and for r] < 3.8 x 10~^°, this corresponds to < 0.08 eV. 

The subsequent non-linear growth in Sp/p was discussed in these sessions in some detail at 
this meeting. In particular, there was a considerable discussion of the effects of the non-linear 
regime on the power spectrum and the appearance of non-Gaussian features such as skewness 
and kurtosis [^. Colombi discussed the 2- and 3-point correlation functions. Numerous 



simulations were presented to reflect the hydrodynamics of galaxy formation [0] and the origin 
of the large scale bias . 



Because we are considering adiabatic perturbations, there will be anisotropies produced 
in the microwave background radiation on the order of 6T/T ~ 6. The value of S, the amplitude 



of the density fluctuations at horizon crossing, has now been determined by COBE 6 = 
(5.7 ± 0.4) X 10"^. Without the existence of dark matter, 6p/p in baryons could then achieve 
a maximum value of only Sp/p ~ AxS{T„i/To) ^ 2 x 10^'^ Ax, where To = 2.35 x 10~^ eV is the 
present temperature of the microwave background and ~ 1 — 10 is a scale dependent growth 
factor. The overall growth in Sp/p is too small to argue that growth has entered a nonlinear 
regime needed to explain the large value (10^) of Sp/p in galaxies. 

Dark matter easily remedies this dilemma in the following way. The transition to mat- 
ter dominance is determined by setting equal to each other the energy densities in radiation 
(photons and any massless neutrinos) and matter (baryons and any dark matter). While if we 
suppose that there exists dark matter with an abundance = n^/n^ (the ratio of the number 
density of x's to photons) then 

Tm = 0.22m^Y^ (10) 

Since we can write m^Y^/GeV = f2^/i^/(4 x 10"^), we have T^/To = 2.4 x lO'^Q^h'^ which is 
to be compared with 350 in the case of baryons alone. The baryons, although still bound 
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to the radiation until decoupling, now see deep potential wells formed by the dark matter 
perturbations to fall into and are no longer required to grow at the rate 6p/p (x R. 

With regard to dark matter and galaxy formation, all forms of dark matter are not equal. 



They can be distinguished by their relative temperature at Particles which are still 

largely relativistic at (like neutrinos or other particles with < 100 eV) have the property 
|]28[| that (due to free streaming) they erase perturbations out to very large scales given by the 
Jean's mass 

Mj = 3x 10^^^^ (11) 
"^jy (ev) 

Thus, very large scale structures form first and galaxies are expected to fragment out later. 
Particles with this property are termed hot dark matter particles. Cold particles (m^ > 1 
MeV) have the opposite behavior. Small scale structure forms first aggregating to form larger 
structures later. Neither candidate is completely satisfactory when the resulting structure is 
compared to the observations. For more details, I refer the reader to reviews in refs. [|I|. 

The most promising possibility we have for unscrambling the various possible scenarios 
for structure formation is the careful analysis of the observed power spectrum. Rewriting the 
density contrast in Fourier space, 

S(k) oc / d'^x^{x)e'^-'' (12) 
J p 

the power spectrum is just 

Pik) = {m\') (13) 

which is often written in term of a transfer function and a power law, P ~ T{k)k'^\ (The flat 
spectrum produced by inflation has n = 1.) The data contributing to P{k) include observations 
of galaxy distributions, velocity distributions and of course the CMB |2D|. However to make a 
comparison with our theoretical expectations, we rely on n-body simulations and a deconvolu- 
tion of the theory from the observations. Overall, there is actually a considerable amount of 
concordance with our expectations. The velocity distributions indicate that 0.3 < < 1 and 
the power spectrum is roughly consistent with an f2 = 1, and n = 1, cold dark matter Universe. 



Future surveys |30, 31] will however, most certainly dramatically improve our understand- 
ing of the detailed features of the power spectrum and their theoretical interpretations. We 
should in principle be able to distinguish between a mixed dark matter and a cold dark matter 
f2 = 1 Universe, whether or not Q < 1, with CDM, the presence of a cosmological constant, 
or a tilt in the spectrum. Zurek [^, stressed the importance of numerical simulations in this 



context. Future probes of the small scale anisotropy |33[ should in addition be able to 
determine the values of not only Q, but also Qb and ho to a high degree of accuracy through a 
careful analysis of the Doppler peak in the power spectrum. 

Accepting the dark matter hypothesis, the first choice for a candidate should be something 
we know to exist, baryons. Though baryonic dark matter can not be the whole story if f2 = 1, 
the identity of the dark matter in galactic halos, which appear to contribute at the level of 
Q ~ 0.05, remains an important question needing to be resolved. A baryon density of this 
magnitude is not excluded by nucleosynthesis. Indeed we know some of the baryons are dark 
since Q < 0.01 in the disk of the galaxy. 



It is quite difficult, however, to hide large amounts of baryonic matter [34]. Sites for halo 



baryons that have been discussed include snowballs, which tend to sublimate, cold hydrogen 
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gas, which must be supported against collapse, and hot gas, which can be excluded by the 
X-ray background. Stellar objects (collectively called MACHOs for macroscopic compact halo 
objects) must either be so small ( M < 0.08 Mq) so as not to have begun nuclear burning or 
so massive so as to have undergone total gravitational collapse without the ejection of heavy 
elements. Intermediate mass stars are generally quite problematic because either they are 
expected to still reside on the main-sequence today and hence would be visible, or they would 
have produced an excess of heavy elements. 

On the other hand, MACHOs are a candidate which are testable by the gravitational 
microlensing of stars in a neighboring galaxy such as the LMC [^. By observing millions 
of stars and examining their intensity as a function of time, it may be possible to determine 
the presence of dark objects in our halo. It is expected that durin g a lensing event, a star's 
intensity will rise in an achromatic fashion over a period 5t ~ 3 \J M/.OOIM© days. Indeed, 
microlensing candidates have been found For low mass objects, those with M < O.IMq, 



it appears however that the halo fraction of Machos is small, ~ 0.191 j^g El- There have been 
some recent reports of events with longer duration |^ leading to the speculation of a white 



dwarf population in the halo. Though it is too early to determine the implications of these 
observations, they are very encouraging in that perhaps this issue can and will be decided. 

The degree to which baryons can contribute to dark matter depends ultimately on the 
overall baryon contribution to Q which is constrained by nucleosynthesis. Because of its impor- 
tance regarding the issue of dark matter and in particular non-baryonic dark matter, I want to 
review the status of big bang nucleosynthesis. 

Conditions for the synthesis of the light elements were attained in the early Universe at 
temperatures T < 1 MeV, corresponding to an age of about 1 second. Given a single input 
parameter, the baryon-to-photon ratio, 1], the theory is capable of predicting the primordial 
abundances of D/H, '^He/H, ^He/H and ^Li/H. The comparison of these predictions with the 
observational determination of the abundances of the light elements not only tests the theory 
but also fixes the value of r]. 

The dominant product of big bang nucleosynthesis is ^He resulting in an abundance of 
close to 25 % by mass. In the standard model, the "^He mass fraction depends only weakly on 
rj. When we go beyond the standard model, the ^He abundance is very sensitive to changes 
in the expansion rate which can be related to the effective number of neutrino flavors. Lesser 
amounts of the other light elements are produced: D and '^He at the level of about 10~^ by 
number, and ^Li at the level of 10~^° by number. 

There is now a good collection of abundance information on the ^He mass fraction, Y, 



0/H, and N/H in over 50 extragalactic HII (ionized hydrogen) regions The obser- 

vation of the heavy elements is important as the helium mass fraction observed in these HII 
regions has been augmented by some stellar processing, the degree to which is given by the 
oxygen and nitrogen abundances. In an extensive study based on the data in []39| , |40[| , it was 
found that the data is well represented by a linear correlation for Y vs. 0/H and Y vs. 



N/H. It is then expected that the primordial abundance of ^He can be determined from the 
intercept of that relation. The overall result of that analysis indicated a primordial mass frac- 
tion, Yp = 0.232 ± 0.003. In ||4^, the stability of this fit was verified by a Monte-Carlo analysis 



showing that the fits were not overly sensitive to any particular HII region. In addition, the 
data from was also included, yielding a "^He mass fraction 



yp = 0.234 ±0.003 ±0.005 (14) 
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The second uncertainty is an estimate of the systematic uncertainty in the abundance determi- 
nation. 

The ''Li abundance is also reasonably well known. In old, hot, population-II stars, ^Li 
is found to have a very nearly uniform abundance. For stars with a surface temperature 
T > 5500 K and a metallicity less than about l/20th solar (so that effects such as stellar 
convection may not be important), the abundances show little or no dispersion beyond that 
which is consistent with the errors of individual measurements. Indeed, as detailed in |Q , much 



of the work concerning ^Li has to do with the presence or absence of dispersion and whether 
or not there is in fact some tiny slope to a [Li] = log "^Li/H + 12 vs. T or [Li] vs. [Fe/H] 
relationship. There is ''Li data from nearly 100 halo stars, from a variety of sources. I will 



use the value given in [45 as the best estimate for the mean "^Li abundance and its statistical 
uncertainty in halo stars 

Li/H = (1.6 ±0.1^°:™)x 10-1° (15) 

The first error is statistical, and the second is a systematic uncertainty that covers the range 
of abundances derived by various methods. The third set of errors in Eq. ([15|) accounts for 
the possibility that as much as half of the primordial ''Li has been destroyed in stars, and that 
as much as 30% of the observed ^Li may have been produced in cosmic ray collisions rather 
than in the Big Bang. Observations of ^Li, Be, and B help constrain the degree to which these 
effects play a role [HH]. For ''Li, the uncertainties are clearly dominated by systematic effects. 



Turning to D/H, we have three basic types of abundance information: 1) ISM data, 
2) solar system information, and perhaps 3) a primordial abundance from quasar absorption 
systems. The best measurement for ISM D/H is [^ 



(D/H)isM = 1.60 ± O.OQ^Ho X 10"^ (16) 

This value may not be universal (or galactic as the case may be) since there may be some real 
dispersion of D/H in the ISM The solar abundance of D/H is inferred from two distinct 
measurements of ^He. The solar wind measurements of ^He as well as the low temperature 
components of step-wise heating measurements of ^He in meteorites yield the presolar (D + 
^He) /H ratio, as D was efficiently burned to ^He in the Sun's pre-main-sequence phase. These 
measurements indicate that [H^, 



D+ ^He 



(4.1±0.6±1.4) X 10~^ (17) 



V H 

The high temperature components in meteorites are believed to yield the true solar '^He/H ratio 

/ ^He^ 





3l 



of Eg, 50 



H 



I = (1.5 ±0.2 ±0.3) X 10-^ 



The difference between these two abundances reveals the presolar D/H ratio, giving. 



(D/H)q ^ (2.6±0.6±1.4) X 10"^ (19) 

Finally, there have been several recent reported measurements of D /H is high redshift quasar 
absorption systems. Such measurements are in principle capable of determining the primordial 
value for D/H and hence t], because of the strong and monotonic dependence of D/H on rj. 
However, at present, detections of D/H using quasar absorption systems indicate both a high 
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and low value of D/H. As such, we caution that these values may not turn out to represent 
the true primordial value. The first of these measurements indicated a rather high D/H 
ratio, D/H ^ 1.9 - 2.5 xlO~^. A recent re-observation of the high D/H absorption system has 
been resolved into two components, both yielding high values with an average value of D /H = 
1.9to;3 ^ 10~^ [0 as well as an additional system with a similar high value |5^. Other high 
D/H ratios were reported in However, there are reported low values of D/H in other such 
systems with values D/H ~ 2.5 x 10~^, significantly lower than the ones quoted above. It 
is probably premature to use this value as the primordial D/H abundance in an analysis of big 
bang nucleosynthesis, but it is certainly encouraging that future observations may soon yield 
a firm value for D/H. It is however important to note that there does seem to be a trend that 
over the history of the Galaxy, the D/H ratio is decreasing, something we expect from galactic 
chemical evolution. Of course the total amount of deuterium astration that has occurred is still 
uncertain, and model dependent. 

There are also several types of ^He measurements. As noted above, meteoritic extractions 
yield a presolar value for ^He/H as given in Eq. ([18|) . In addition, there are several ISM 
measurements of ^He in galactic HII regions which also show a wide dispersion 



H 



~ 1 - 5 X 10"^ 



(20) 



HII 



There is also a recent ISM measurement of "^He 15711 with 




2.lti X 10" 



(21) 



ISM 



Finally there are observations of He in planetary nebulae [58| which show a very high He 



abundance of ^He/H ~ 10~^. 

Each of the light element isotopes can be made consistent with theory for a specific range 
in 1]. Overall consistency of course requires that the range in t] agree among all four light 
elements. ^He (together with D) has stood out in its importance for BBN, because it provided 
a (relatively large) lower limit for the baryon-to-photon ratio rjiQ > 2.8. This limit for 
a long time was seen to be essential because it provided the only means for bounding t] from 
below and in effect allows one to set an upper limit on the number of neutrino flavors 



Njy, as well as other constraints on particle physics properties. That is, the upper bound to 
is strongly dependent on the lower bound to rj. This is easy to see: for lower rj, the ^He 
abundance drops, allowing for a larger N^, which would raise the ^He abundance. However, 
for 1] < 4 X 10"^^, corresponding to Qh"^ ~ .001 — .002, which is not too different from galactic 
mass densities, there is no bound whatsoever on N,, Of course, with the improved data 
on ^Li, we do have lower bounds on r] which exceed 10~^°. 

In it was argued that since stars (even massive stars) do not destroy '^He in its 
entirety, we can obtain a bound on t] from an upper bound to the solar D and ^He abundances. 
One can in fact limit ||5^, |6^ the sum of primordial D and ^He by applying the expression below 
at t = 

'D + 3He\ /D\ 1 /3He\ , , 



H 



In (^2]) , is the fraction of a star's initial D and ^He which survives as ^He. For > 0.25 as 
suggested by stellar models, and using the solar data on D/H and '^He/H, one finds r/ig > 2.8. 
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The limit t]iq > 2.8 derived using (^) is really a one shot approximation. Namely, it 
is assumed that material passes through a star no more than once. To determine whether 
or not the solar (and present) values of D/H and '^He/H can be matched it is necessary to 
consider models of galactic chemical evolution |Q. In the absence of stellar ^He production, 
particularly by low mass stars, it was shown that there are indeed suitable choices for a star 
formation rate and an initial mass function to: 1) match the D/H evolution from a primordial 
value (D/H)p = 7.5 x 10~^, corresponding to r^io = 3, through the solar and ISM abundances, 
while 2) at the same time keeping the ^He/H evolution relatively flat so as not to overproduce 
■^He at the solar and present epochs. This was achieved for > 0.3. Even for ~ 0.7, the 
present ^He/H could be matched, though the solar value was found to be a factor of 2 too high. 
For (D/H)p ~ 2 X 10~^, corresponding to rjiQ ~ 1.7, models could be found which destroy D 
sufficiently; however, overproduction of '^He occurred unless g^ was tuned down to about 0.1. 

In the context of models of galactic chemical evolution, there is, however, little justification 
a priori for neglecting the production of ^He in low mass stars. Indeed, stellar models predict 
that considerable amounts of ^He are produced in stars between 1 and 3 M0. For M < 8Mq, 
Iben and Truran calculate 

(^He/H)j = 1.8 X 10-M ^ j +0.7[(D+ =^He)/H]. (23) 

so that at r]io = 3, and ((D + 3He)/H)i = 9 x 10~^ 5-3(11^10) = 2.7! It should be emphasized 
that this prediction is in fact consistent with the observation of high ^He/H in planetary nebulae 
|58|. 

Generally, implementation of the ^He yield in Eq. ( p3D in chemical evolution models 



leads to an overproduction of ^He/H particularly at the solar epoch [^, In Figure 1, the 
evolution of D/H and ^He/H is shown as a function of time from p9| , |6^ . The solid curves 
show the evolution in a simple model of galactic chemical evolution with a star formation rate 
proportional to the gas density and a power law IMF (see [|^) for details). The model was 
chosen to fit the observed deuterium abundances. However, as one can plainly see, '^He is grossly 
overproduced (the deuterium data is represented by squares and '^He by circles). Depending 
on the particular model chosen, it may be possible to come close to at least the upper end of 
the range of the ^He/H observed in galactic HII regions [^, although the solar value is missed 
by many standard deviations. 

The overproduction of ^He relative to the solar meteoritic value seems to be a generic 
feature of chemical evolution models when ^He production in low mass stars is included. In 
|9|, a more extreme model of galactic chemical evolution was tested. There, it was assumed 



that the initial mass function was time dependent in such a way so as to favor massive stars early 
on (during the first two Gyr of the galaxy). Massive stars are preferential from the point of view 
of destroying ^He. However, massive stars are also proficient producers of heavy elements and 
in order to keep the metallicity of the disk down to acceptable levels, supernovae driven outflow 
was also included. The degree of outflow was limited roughly by the observed metallicity in 
the intergalactic gas in clusters of galaxies. One further assumption was necessary; we allowed 
the massive stars to lose their '^He depleted hydrogen envelopes prior to explosion. Thus 
only the heavier elements were expulsed from the galaxy. With all of these (semi-defensible) 
assumptions, ^He was still overproduced at the solar epoch, as shown by the dashed curve in 
Figure 1. Though there certainly is an improvement in the evolution of ^He without reducing 
the yields of low mass stars, it is hard to envision much further reduction in the solar '^He 
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predicted by these models. The only conclusion that we can make at this point is that there is 
clearly something wrong with our understanding of ^He in terms of either chemical evolution, 
stellar evolution or perhaps even the observational data. 

Given the magnitude of the problems concerning ^He, it would seem unwise to make any 
strong conclusion regarding big bang nucleosynthesis which is based on ^He. Perhaps as well 
some caution is deserved with regard to the recent D/H measurements, although if the present 
trend continues and is verified in several different quasar absorption systems, then D/H will 
certainly become our best measure for the baryon-to-photon ratio. Given the current situation 
however, it makes sense to take a step back and perform an analysis of big bang nucleosynthesis 
in terms of the element isotopes that are best understood, namely ^He and ^Li. 

Monte Carlo techniques are proving to be a useful form of analysis for big bang nucle- 
osynthesis [|68| , |69| . In [|I8|, we performed just such an analysis using only ^He and ''Li. It 



should be noted that in principle, two elements should be sufficient for constraining the one 
parameter (77) theory of BBN. We begin by establishing likelihood functions for the theory and 
observations. For example, for ^He, the theoretical likelihood function takes the form 

Lbb^{Y, Fbbn) = e-(^-^BBN(.))^/2.? (24) 

where Ybbn{v) is the central value for the ^He mass fraction produced in the big bang as 
predicted by the theory at a given value of 1], and o"i is the uncertainty in that value derived 
from the Monte Carlo calculations and is a measure of the theoretical uncertainty in 



the big bang calculation. Similarly one can write down an expression for the observational 
likelihood function. In this case we have two sources of errors, as discussed above, a statistical 
uncertainty, a2 and a systematic uncertainty, cXsys. Here, I will assume that the systematic error 



is described by a top hat distribution [^, |70|. The convolution of the top hat distribution and 



the Gaussian (to describe the statistical errors in the observations) results in the difference of 
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0.00 2.00 4.00 6.00 8.00 10.00 

Figure 2: Likelihood distribution for each of ^He and ^Li, shown as a function of rj. The 
one-peak structure of the ^He curve corresponds to its monotonic increase with 77, while the 
two-peaks for ^Li arise from its passing through a minimum. 



two error functions 



Y-Yo + a,yA r Y-Yo 



yo) ^ erf [^1^) - erf [^tr") (25) 

where in this case, Yq is the observed (or observationally determined) value for the ^He mass 
fraction. (Had I used a Gaussian to describe the systematic uncertainty, the convolution of 
two Gaussians leads to a Gaussian, and the likelihood function (UBi) would have taken a form 



similar to that in (p^). 

A total likelihood function for each value of 7710 is derived by convolving the theoretical 
and observational distributions, which for ^He is given by 

L'^^'toUr]) = J dYLBBN {Y, Fbbn (v)) Lo{Y, Fq) (26) 

An analogous calculation is performed ||T^ for ^Li. The resulting likelihood functions from the 
observed abundances given in Eqs. (0) and ( plSj) is shown in Figure 2. As one can see there is 
very good agreement between ^He and ^Li in the vicinity of rjio ~ 1.8. 

The combined likelihood, for fitting both elements simultaneously, is given by the product 
of the two functions in Figure |^ and is shown in Figure ^. From Figure it is clear that 
^He overlaps the lower (in rj) ^Li peak, and so one expects that there will be concordance 
in an allowed range of r] given by the overlap region. This is what one finds in Figure ^ 
which does show concordance and gives a preferred value for 7], t]iq = 1.8ll\ corresponding to 

n^h^ = meliZv 

Thus, we can conclude that the abundances of ^He and ^Li are consistent, and select 
an rjiQ range which overlaps with (at the 95% CL) the longstanding favorite range around 
rjiQ = 3. Furthermore, by finding concordance using only *^He and ''Li, we deduce that if there 
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Figure 3: Combined likelihood for simultaneously fitting ^He and ^Li, as a function of rj. 



is problem with BBN, it must arise from D and ^He and is thus tied to chemical evolution or the 
stellar evolution of '^He. The most model-independent conclusion is that standard BBN with 
Ni, = 3 is not in jeopardy, but there may be problems with our detailed understanding of D and 
particularly ^He chemical evolution. It is interesting to note that the central (and strongly) 
peaked value of rjio determined from the combined ^He and ''Li likelihoods is at rjio = 1.8. 
The corresponding value of D/H is 1.8 xlO~^, very close to the high value of D/H in quasar 



absorbers |51, 52, 54|. Since D and ^He are monotonic functions of rj, a prediction for t], based 
on "^He and ^Li, can be turned into a prediction for D and ^He. The corresponding 95% CL 
ranges are D/H = (5.5 - 27) x 10"^ and and ^He/H = (1.4 - 2.7) x 10"^ 

If we did have full confidence in the measured value of D/H in quasar absorption systems, 
then we could perform the same statistical analysis using "^He, ^Li, and D. To include D/H, 
one would proceed in much the same way as with the other two light elements. We compute 
likelihood functions for the BBN predictions as in Eq. (^) and the likelihood function for the 
observations using D/H = (1.9 ± 0.4) x 10~^. We are using only the high value of D/H here. 
These are then convolved as in Eq. (^6]). In figure 4, the resulting normalized distribution, 
^uitaiiv) is super-imposed on distributions appearing in figure 2. It is indeed startling how the 
three peaks, for D, ^He and ^Li are literally on top of each other. In figure 5, the combined 
distribution is shown. We now have a very clean distribution and prediction for r], r]iQ = 1.751;^ 
corresponding to fi^^^ = •006t;Q0Q4, with the peak of the distribution at rjiQ = 1.75. The 
absence of any overlap with the high-?] peak of the ^Li distribution has considerably lowered 
the upper limit to rj. Overall, the concordance limits in this case are dominated by the deuterium 
likelihood function. 

To summarize on the subject of big bang nucleosynthesis, I would assert that one can 
conclude that the present data on the abundances of the light element isotopes are consistent 
with the standard model of big bang nucleosynthesis. Using the the isotopes with the best 
data, ''He and ^Li, it is possible to constrain the theory and obtain a best value for the baryon- 
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Figure 4: As in Figure 2, with the addition of the hkehhood distribution for D/H. 
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Figure 5: Combined hkehhood for simultaneously fitting ^He and ^Li, and D as a function of 
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to-photon ratio of rjio = 1.8, a corresponding value fi^ = .0065 and 



1.4 < r]w < 3.8 95%CL 
.005 < nBh^ < .014 95%CL (27) 

For 0.4 < h < 1, we have a range .005 < Qb < -09. This is a rather low value for the baryon 



density and would suggest that much of the galactic dark matter is non-baryonic |71|. These 
predictions are in addition consistent with recent observations of D/H in quasar absorption 
systems which show a high value. Difficulty remains however, in matching the solar ^He abun- 
dance, suggesting a problem with our current understanding of galactic chemical evolution or 
the stellar evolution of low mass stars as they pertain to ^He. 

If we now take as given that non-baryonic dark matter is required, we are faced with the 
problem of its identity. Light neutrinos (m < 30eV^) are a long-time standard when it comes 



to non-baryonic dark matter |72|. Light neutrinos produce structure on large scales and the 
natural (minimal) scale for structure clustering is given in Eq. (pA]). Hence neutrinos offer the 
natural possibility for large scale structures ||73|, WM including filaments and voids. It seemed. 



however, that neutrinos were ruled out because they tend to produce too much large scale 



structure |7^. Because the smallest non- linear structures have mass scale Mj and the typical 
galactic mass scale is ~ IO^^Mq, galaxies must fragment out of the larger pancake-like objects. 



The problem is that in such a scenario, galaxies form late |7^, |76| {z < 1) whereas quasars and 



galaxies are seen out to redshifts z > A. Recently, neutrinos are seeing somewhat of a revival 
in popularity in mixed dark matter models. 

In the standard model, the absence of a right-handed neutrino state precludes the existence 
of a neutrino mass. By adding a right-handed state ur, it is possible to generate a Dirac mass 
for the neutrino, = h,^v/\/2, as for the charged lepton masses, where h^, is the neutrino 
Yukawa coupling constant, and v is the Higgs expectation value. It is also possible to generate 
a Majorana mass for the neutrino when in addition to the Dirac mass term, m^ujii'L, a term 
Mv^vji is included. In what is known as the see-saw mechanism, the two mass eigenstates are 
given by niy-^ ~ m1/M which is very light, and ~ M which is heavy. The state vi is our 
hot dark matter candidate as z/2 is in general not stable. 

The cosmological constraint on the mass of a light neutrino is derived from the overall 
mass density of the Universe. In general, the mass density of a light particle x can be expressed 
as 

= m^Y^n^ <pc = 1-06 x lO-^/io^GeV/cm^ (28) 
where Y^ = n^/uy is the density of x's relative to the density of photons, for ilho^ < 1. For 



neutrinos Yy = 3/11, and one finds |]77| 



Y^fyrriy < dSeViflK') (29) 
u 2 

where the sum runs over neutrino fiavors. All particles with abundances Y similar to neutrinos 
will have a mass limit given in Eq. (p9|). 

It was possible that neutrinos (though not any of the known fiavors) could have had large 



> 1 MeV. In that case their abundance Y is controlled by z/, z/ annihilations [^, 
for example, uu —>■ ff via Z exchange. When the annihilations freeze-out (the annihilation 
rate becomes slower than the expansion rate of the Universe), Y becomes fixed. Roughly, 
Y ~ (mcr^)~^ and p ~ cr^^^ where is the annihilation cross-section. For neutrinos, we 
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^'jmz^ so that ~ ^jva^ and we can derive a lower bound |7^, ^ on the 



expect a A ~ 

neutrino mass, niy > 3 — 7 GeV, depending on whether it is a Dirac or Majorana neutrino. 
Indeed, any particle with roughly a weak scale cross-sections will tend to give an interesting 
value of ^h? ' 



Due primarily to the limits from LEP |8^, the heavy massive neutrino has become simply 
an example and is no longer a dark matter candidate. LEP excludes neutrinos (with standard 
weak interactions) with masses niy < 40 GeV. Lab constraints for Dirac neutrinos are available 
|83[| , excluding neutrinos with masses between 10 GeV and 4.7 TeV. This is significant, since 
it precludes the possibility of neutrino dark matter based on an asymmetry between u and u 
Majorana neutrinos are excluded as dark matter since n^ho^ < 0.001 for > 40 GeV 



and are thus cosmologically uninteresting. 

Supersymmetric theories introduce several possible candidates. If R-parity, which distin- 
guishes between "normal" matter and the supersymmetric partners and can be defined in terms 
of baryon, lepton and spin as R = (—1)3^+^+25'^ jg unbroken, there is at least one supersym- 
metric particle which must be stable. I will assume R-parity conservation, which is common 
in the MSSM. R-parity is generally assumed in order to justify the absence of superpotential 
terms can be responsible for rampid proton decay. The stable particle (usually called the LSP) 
is most probably some linear combination of the only R = —1 neutral fermions, the neutralinos 
85[| : the wino W^, the partner of the 3rd component of the SU{2)l gauge boson; the bino. 



B, the partner of the U{1)y gauge boson; and the two neutral Higgsinos, Hi and H2. Gluinos 
are expected to be heavier -rrig = sin^ 6wM2, where M2 is the supersymmetry breaking 
SU(2) gaugino mass- and they do not mix with the other states. The sneutrino 
possibility but has been excluded as a dark matter candidate by direct searches 



is also a 
indirect 

8% and accelerator |^2|] searches. For more on the motivations for supersymmetry and the 



supersymmetric parameter space, see the contribution of Jungman pT|. 

The identity of the LSP is effectively determined by three parameter in the MSSM, the 
gaugino mass, M2, the Higgs mixing mass n, and the ratio of the Higgs vacuum expectation 
values, tan/3. In Figure 6 regions in the M2, fi plane with tan/3 = 2 are shown in which 
the LSP is one of several nearly pure states, the photino, 7, the U(l) gaugino, B, a symmetric 
combination of the Higgsinos, H{i2) = ■^(-^i+-^2), or the Higgsino S = Hi cos[3+H2 sin j3. The 
dashed lines show the LSP mass contours. The cross hatched regions correspond to parameters 
giving a chargino {W^,H^) state with mass rriy. < AbGeV and as such are excluded by LEPI 



This constraint has been extended by LEPI. 5, [^T], |9^ and is shown by the light shaded region 
and corresponds to regions where the chargino mass is < 67 GeV. The dark shaded region 
corresponds to a limit on M2 from the limit on the gluino mass rrig < 70 GeV or M2 < 22 
GeV. Notice that the parameter space is dominated by the B or H12 pure states and that the 
photino (most often quoted as the LSP) only occupies a small fraction of the parameter space, 
as does the Higgsino combination S^. 

As described in |^T|, the relic abundance of LSP's is determined by solving the Boltzmann 
equation for the LSP number density in an expanding Universe. The technique ||80|| used is 
similar to that for computing the relic abundance of massive neutrinos ||78||. For binos, as was 
the case for photinos it is possible to adjust the sfermion masses to obtain closure 
density. Adjusting the sfermion mixing parameters or CP violating phases |^ allows even 
greater freedom. In Figure 7 [^, the relic abundance {Qh"^) is shown in the M2 — /i plane with 
tan/3 = 2, the Higgs pseudoscalar mass mo = 50 GeV, mt = 100 GeV, and fnj = 200 GeV. 
Clearly the MSSM offers sufficient room to solve the dark matter problem. Similar results have 
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Figure 6: The M2-H plane in the MSSM for tan/5 = 2. 



been found by other groups |100 



In Figure 7, in the higgsino sector H12 marked off by 
the dashed hne, co-annihilations [ p.01| , |99|] between -ff(i2) and the next lightest neutralino (also 
a Higgsino) were not included. These tend to lower significantly the relic abundance in much 
of this sector. 

Though I have concentrated on the LSP in the MSSM as a cold dark matter candidate, 
there are many other possibilities when one goes beyond the standard model. Axions were 
discussed at length by Jungman ||8T| and Lu [|102|| . A host of other possibihties were discussed 
by Khlopov [|Tol . 



The final subject that I will cover in this introduction/summary is the question of detec- 
tion. Dark matter detection can be separated into two basic methods, direct ||104|| and indirect 
TUB. Direct detection relies on the ability to detect the elastic scattering of a dark matter 



candidate off a nucleus in a detector. The experimental signatures for direct detection were 
covered by Cabrera | 1U6 | and several individual experiments were described | 1U7 . 

The detection rate, will depend on the density of dark matter in the solar neighborhood, 
p ~ 0.3 GeV/cm^, the velocity, v ~ 300 km/s, and the elastic cross section, a. Spin indepen- 
dent interactions are the most promising for detection. Dirac neutrinos have spin-independent 
interactions, but as noted above, these have already been excluded as dark matter by direct 
detection experiments [p3i. In the MSSM, it is possible for the LSP to also have spin indepen- 



dent interactions which are mediated by Higgs exchange. These scatterings are only important 
when the LSP is a mixed (gaugino/Higgsino) state as in the central regions of Figures 6 and 7. 
Generally, these regions have low values of Qh"^ (since the annihilation cross sections are also 
enhanced) and the parameter space in which the elastic cross section and relic density are large 
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Figure 7: Relic density contours (fi/i^) in the M2 - /x plane. 



is rather limited. Furthermore, a significant detection rate in this case relies on a low mass for 
the Higgs scalar [|Tn|, M . 

More typical of the SUSY parameter space is a LSP with spin dependent interactions. 
Elastic scatterings are primarily spin dependent whenever the LSP is mostly either gaugino or 
Higgsino. For Higgsino dark matter, Higgsinos with scatterings mediated by Z° avoid the -f^(i2) 
regions of Figures 6 and 7, and as such are now largely excluded (the region does grow at 
low tan/3 |8^. Higgsino scatterings mediated by sfermion exchange depend on couplings 
proportional to the light quark masses and will have cross sections which are suppressed by 
(mp/m^y)^, where rrip is the proton mass. These rates are generally very low ||109|| . Binos, on 
the other hand, will have elastic cross sections which go as m^/mj^, where m is the reduced 
mass of the bino and nucleus. These rates are typically higher (reaching up to almost 0.1 events 
per kg-day [|10|, [Tig, [TTT |. 



Indirect methods also offer the possibility for the detection of dark matter. Three methods 
for indirect detection were discussed [ [105|| . 1) 7-rays from dark matter annihilations in the 
galactic halo are a possible signature ||112|| . In the case of the MSSM, unless the mass of the LSP 
is larger than mw, the rates are probably too small to be detectable over background ||105|| . 2) 
Dark matter will be trapped gradually in the sun, and annihilations within the sun will produce 
high energy neutrinos which may be detected ; similarly, annihilations within the earth may 
provide a detectable neutrino signal ||1 14|] . Edsjo [|115|| discussed possibilities for determining 
the mass the dark matter candidate from the angular distribution of neutrinos. This method 
hold considerable promise, as there will be a number of very large neutrino detectors coming 
on line in the future. Finally, 3) there is the possibility that halo annihilations into positrons 
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and antiprotons in sufficient numbers to distinguish them from cosmic-ray backgrounds ||112 
riel ITTTll. 
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